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The  proposed  research  was  to  apply  the  principles 
of  computer  aid  tomography  to  acousto-optic  data  for  the 
purpose  of  investigating  sound  radiation  and  scattering. 
The  principal  investigator  had  previously  developed 
the  theory  and  logic  but  lacked  the  facilities,  namely 
for  computer  acquisition,  mass  storage,  display  of  the 
large  amounts  of  information  involved.  The  objective 
of  this  contract  was  to  allow  the  principal  investigator 
and  colleagues  at  the  Navy  Research  Laboratory, 

Washington,  DC  to  implement  the  necessary  hardware  and 
computer  algorithms  during  the  semester  of  1979.  The 
Physical  Acoustics  Branch  of  the  Naval  Research  Laboratory 
had  both  the  acousto-optic  facilities  and  the  on-line  com¬ 
puter  system.  Upon  submission  of  the  proposal,  the 
principal  investigator  sent  to  the  scientific  personnel  at 
NRL,  the  requirements  for  the  manipulation  of  the  acoustic 
radiating  system.  The  personnel  there  designed  a 
manipulator  and  ordered  the  required  parts  to  build 
the  device. 

The  plan  of  action  was  for  the  principal  investigator 
to  become  familar  with  the  computer  system  while  the 
manipulator  was  being  built.  The  parts  required  for  the 
manipulator  were  not  received  during  the  summer  of  1979. 

Consequently,  the  Principal  Investigator  devoted 
himself  to  theoretical  and  computational  problems  of 
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experiments  underway  at  URL  pertaining  to  acou st o -opt ic s 
and  radiation  of  sound  fields. 

ACOUSTO-OPTICS 

In  as  much  as  tomography  may  be  considered  as  the 
inversion  of  a  scattering  process,  the  Principal  Investi¬ 
gator  developed  an  algorithm  to  calculate  the  direct 
process  of  scattering  of  light  by  ultrasound.  The 
algorithm  is  novel  in  that,  one  needs  only  to  specify  the 
nature  of  radiating  acoustic  transducer  and  many  of  the 
observed  acousto-optic  phenomena  can  be  predicted  and 
graphically  displayed.  Attachment  A  and  B  are  publications 
resulting  from  this  activity. 

FOCUSED  ULTRASONIC  FIELDS 

An  ongoing  activity  on  the  Physical  Acoustic  Branch 
of  NRL  was  experimental  exploration  of  focused  ultrasonic 
fields.  The  Principal  Investigator  had  recently  super¬ 
vised  a  PhD  student  on  the  theory  of  such  fields.  A 
joint  experimental  and  computational  effort  was  made 
to  validate  the  theory.  Attachment  C  is  a  publication 
which  shows  the  success  of  the  activity. 

ACOUSTO-OPTIC  TOMOGRAPHY 

The  contract  was  extended  to  May  1980  to  allow  the 
Principal  Investigator  to  publish  the  material  and 


continue  his  theoretical  investigation  at  the  University 
of  Houston.  During  this  period  he  developed  a  new 
theory  of  data  acquisition  for  acousto-optic  tomography, 
trained  two  undergraduates  and  submitted  a  new  proposal 
to  the  Office  of  Naval  Research  in  addition  to  publishing. 
The  new  proposal  was  funded  by  ONR  and  the  two  students 
and  the  Principal  Investigator  returned  to  NRL  for  the 
Summer  of  1980  to  continue  the  research. 
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A  Numerical  Procedure  for  Calculating  the 
Integrated  Acoustooptic  Effect  £ 

BILL  D.  COOK,  senior  member,  ieee,  EDUARDO  CAVANAGH,  and  HENRY  D.  DARDY,  member,  ieee 


Abstract- The  integrated  optical  effect  of  light  passing  through  an 
ultrasonic  field  can  be  numerically  calculated  from  the  velocity  dis¬ 
tribution  assigned  to  a  planar  transducer.  The  theory  is  developed  for 
linear  acoustics,  nonabsorbing  media,  and  normal  incidence  of  light. 
As  a  consequence  of  the  application  of  the  Fourier  projeclion  slice 
theorem,  an  algorithm  based  on  digital  Fourier  transforms  calculates 
the  integrated  optical  effect  without  integral  ion.  Demonstrations  of 
the  influence  of  transducer  geometry  are  shown  for  the  integrated 
optical  effect  and  schliercn  patterns  based  on  this  effect. 


Introduction 


THE  OPTICAL  methods  for  estimating  field  quantities, 
such  as  pressure,  of  ultrasonic  fields  required  a  model  of 
the  sound  field.  Often  these  were  simplistic  (collimated  plane 
waves)  but  very  useful.  In  reality  the  light  usually  passes 
through  the  near-field  region  of  the  ultrasonic  transducer.  The 
multitude  of  schlieren  photographs  in  the  literature  give  testi¬ 
mony  to  the  complexity  of  the  effect  of  the  near  field.  For 
example.  Fig.  I  shows  a  schlieren  image,  published  by  Oster- 
hammel  (I )  in  1941,  of  the  near  field  of  a  square  transducer. 
The  intricate  pattern  was  explained  by  Osterhammel  as  inter¬ 
ference  of  planes  waves  front  the  face  of  transducer  and  circu¬ 
lar  waves  front  the  edge;  this  theory  does  not  permit  quantita¬ 
tive  calculations  and  is  not  extendable  to  other  geometries. 

Schlieren  patterns  such  as  is  shown  in  Fig.  1  can  be  com¬ 
puted  provided  the  field  parameters  of  the  integrated  optical 
effect  are  known.  At  low  ultrasonic  frequencies,  the  interac¬ 
tion  of  sound  and  light  can  be  described  by  one  parameter, 
the  so-called  Raman-Nath  parameter  v  [2] .  For  a  three- 
dimensional  pressure  field  p(x,y,  z)  etu>t,  with  the  light 
traversing  in  the  y  direction  as  shown  in  Fig.  2,  the  integrated 
optical  effect  is  defined  as 


v  (x,  z) 


2ttk  f 

~JP 


( x,y,z)dy . 


(1) 


Here  X  is  the  optical  wavelength  and  k  is  the  piezooptic  coeffi¬ 
cient.  The  integrated  optical  effect  is  thus  a  two-dimensional 
phasor  field  v  (x,  z)  In  this  paper  a  computation  procedure  is 
presented  for  calculating  v  (x,  z)  when  the  motion  of  the  trans¬ 
ducer  is  specified.  The  procedure  is  straightforward  and  eased 
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Fig.  I .  Schlieren  photograph  of  near  field  of  a  square  transducer  with 
ratio  of  width  to  acoustical  wavelength  of  14.37,  from  Osterhammel 
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Fig.  2.  Diagram  showing  geometry. 


by  utilization  of  the  standard  one-dimensional  Fourier  trans- 
form  which  yields  the  field  v  (x,  z)  at  a  given  z  as  a  one-step 
process. 

Analytic  solutions  and  numerical  procedures  have  been  at¬ 
tempted  on  the  problem  previously.  Ingenito  and  Cook  (3] 
have  given  an  exact  formulation  for  piston-like  motions  of  the 
transducers  with  on-axis  values  given  in  approximation  for  a 
circular,  rectangular,  and  square  piston  on  its  diagonal. 
Maloney,  Meltz.and  Gravel  [4]  have  treated  the  square  trans¬ 
ducer  with  the  Fresnel  approximation.  Their  experimental 
measurements  in  solids  agreed  with  theory  for  both  axial  and 
transverse  profiles.  Crandal  [5]  also  made  computations  and 
measurements  with  emphasis  on  phase  measurements  for 
sound  velocity  measurements.  Haran.Cook,  and  Stewart  (6] 
have  used  the  theoretical  results  for  comparing  optical  calibra¬ 
tion  of  a  circular  transducer  with  a  radiation  pressure  tech¬ 
nique.  Riebold  et  a!  (7] ,  (8J  have  also  made  careful  mea- 
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surements  on  the  axis  values  to  compare  with  the  theory  of 
Igenito  and  Cook  [3) .  Berlinghieri  and  Cook  [9]  have  also 
calculated  the  o-field  for  a  circular  piston  to  test  an  inversion 
algorithm  to  estimate  the  pressure  p  (x,y,  z)  from  measured 
values  of  u(x,  z). 

The  procedure  presented  here  allows  calculations  throughout 
the  near  field  including  the  adjacent  region  to  the  transducer 
where  the  pressure  field  and  integrated  optical  effect  can  vary 
significantly  with  a  small  change  of  position.  The  only  limita¬ 
tion  in  calculating  the  complete  field  appear  to  be  those  con¬ 
nected  with  the  normal  utilization  of  digital  Fourier  trans¬ 
forms. 

Theory 

The  procedure  that  is  to  be  developed  is  based  on  the  appli¬ 
cation  of  the  Fourier  projection  slice  theorem  in  conjunction 
with  the  Fourier  technique  of  solving  the  scalar  wave  equa¬ 
tions.  Briefly  the  Fourier  projection  slice  theorem  states  that 
the  one-dimensional  transform 

W{kx)  =  f  w(x)  exp  (ikxx)  dx  (2) 

eo 

of  the  projection  w  (x),  defined  as 

w(x)=  f  f(x,y)  dy  (3) 

«/ _  t  to 

is  equivalent  to  the  two-dimensional  Fourier  transform 
F(kx,  ky )  with  ky  =  0.  In  other  words  it  states 

W(kx)  =  F(kx,ky*  0).  (4) 

Consequently,  a  knowledge  of  a  slice  of  Fourier  domain 
yields  w(x)  upon  the  inverse  transform.  The  theorem  is  to  be 
applied  to  (1 )  which  is  now  recognized  as  the  form  of  a  pro¬ 
jection  integral.  Thus  the  scheme  is  to  find  the  Fourier  slice 
P(kx,ky  =  0;  z)  corresponding  to  the  local  value  p  (x,y,  z)  of 
the  pressure  produced  by  a  given  normal  particle  velocity 
u(x,y,z  =0)c<u,<  of  the  transducer;  and  then  numerically 
compute  the  inverse  transform. 

A  common  method  of  approaching  scalar  radiation  problems 
is  by  resolving  the  field  parameter,  i.e. ,  pressure,  velocity  po¬ 
tential,  etc.,  into  a  system  of  plane  waves  with  two-dimensional 
Fourier  transforms.  For  example,  let 


as  in  the  case  of  linear  nonabsorbing  acoustic  waves.  Sub¬ 
stituting  (6)  into  (7),  one  can  obtain 

~^(k2-  k2x- k2)P  =  0.  (8) 

A  solution  of  (8)  is 

P(kx,ky'z)  =  P(kx,ky,z0)exp  [-ikz{z  -  z0)J  (9) 

where  kz  =  (k2  -  k\  -  ky)'F,  Thus,  this  approach  permits  cal¬ 
culation  of  the  field  in  the  plane  z  from  a  knowledge  of  the 
field  in  the  plane  z0.  The  quantity  P(kx,  ky\  z)  can  be  inter¬ 
preted  as  the  phasor  description  of  plane  waves  with  w  ave 
vector  ( kx ,  ky,kz).  It  is  noted  that  kz  can  be  real,  which 
means  that  the  wave  propagates,  or  imaginary  which  infers  an 
evanescent  wave. 

Our  intermediate  objective  is  to  find  P(kx,  ky  \  z)  from  the 
velocity  distribution  of  the  transducer.  We  started  with  one 
component  of  the  force  equation 


3  p  duz 

bz'~Plu 

(10a) 

=  -iu>uz  (for  harmonic  waves) 

(10b) 

where  the  z  axis  is  in  the  direction  normal  to  the  transducer 
surface  and  p  is  the  density  of  the  medium .  Substituting  (9) 
for  the  pressure  and  a  similar  expression  for  the  particle  ve¬ 
locity  into  (10)  one  can  find 

P(kx  ,ky  -,z)  =  (k/kz)  pcUz(kx,ky,z)  (11) 

where  c  is  the  velocity  of  sound.  Equation  (1 1)  is  the  plane 
wave  relationship  between  pressure  and  particle  velocity;  the 
term  (k/kz)  in  the  secant  of  the  angle  between  the  direction  of 
propagation  of  the  plane  and  the  z  direction  when  kz  is  real. 

At  this  stage  of  the  development,  we  know  how  to  evaluate 
u(x,z)  from  P(kx,  0;  z)  and  to  calculate/*  (kx,  0;z)at  any 
plane  z  from  P(kx,  0;zo)  which  are  related  to  a  Fourier 
domain  description  of  the  z  component  of  the  particle  ve¬ 
locity  uz{kx,  0;z0).  To  obtain  uz(kx,  0;zo)  from  the  given 
particle  velocity  uz  ( x,y ,  z0),  we  again  apply  the  Fourier  pro¬ 
jection  slice  theorem.  We  take  the  projection  of  uz  (x.y.  z0) 
to  be 

uz  (x,  z0)  =  fu2  (x,  y,  z0)  dy  (12) 


P(kx,ky,z)  =  fh  ,y,z)exp  ['i  (kxx  +  k^y)]  dx  dy 


(5) 


and 


p(x,y,z)  =  (\n)2  JJ  P(kx,ky\z) 


If' 


■  exp  [t(k^x  +  kyy)\  dkx  dky  (6) 

be  a  typical  transform  pair.  Here  p(x,y,  z)  is  a  phasor,  as 
previously  defined,  that  satisfies  the  Helmholtz  equation 

VV  +  *V  =  0  (7) 


where  they  direction  is  chosen  to  be  the  direction  of  light 
propagation.  The  one-dimensional  transform  of  Uz  (x,  z0) 
is  Uz  (kx ,  0;  z0)  our  desired  result. 

Thus  the  procedure  for  calculating  the  field  v(x,  z)  is  as 
follows. 

1 )  For  a  given  direction  of  light  propagation  with  respect  to 
transducer  geometry,  calculate  the  projection  of  velocity 
by  (12). 

2)  Numerically  compute  the  one-dimensiona)  Fourier  trans¬ 
form  to  find  Uz  ( kx ,  0;  z0). 

3)  Divide  by  kz  and  the  other  appropriate  constants  to  ob¬ 
tain  P(kx,  0 ;z0)  as  specified  by  (1 1). 

4)  Calculate  P(kx,  0;  z)  at  desired  plane  z  using  (9). 
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r-ig.  3.  Computed  integrated  optical  effect  (normalized  for  ka  =  10. 
(a)  Circular  transducer,  (b)  Rectangular,  (c)  Square  on  diagonal. 

5)  Compute  the  inverse  one-dimensional  transform  to 
evaluate  v(x,z). 

To  calculate  at  alternate  distance  z  from  the  transducer,  re¬ 
peat  steps  4)  and  S). 


I  Os  o  05  I 


- *■  r/R 

Fig.  4.  Computed  local  variation  of  pressure  for  a  circular  piston  with 
ka  -  10,  from  Meyer  and  Newman  [111.  (Courtesy  of  Academic 
Press). 

Results 

Using  the  procedure  described  above  we  have  calculated  the 
integrated  optical  effect  for  a  circular,  rectangular,  and  square 
transducer  on  its  diagonal  with  ka  =  10.  The  dimension  a  is 
the  radius  of  a  circular  transducer,  half  the  length  of  a  side  of 
the  square  transducer  and  half  the  length  of  the  rectangular 
transducer  perpendicular  to  the  direction  of  light  propagation. 
The  field  variations  of  u(x,  z)  for  these  configurations  are 
shown  in  Fig.  3.  Assuming  the  rectangular  transducer  is  square 
the  normalization  v0  =  4mia/\  has  been  applied  for  all  of  these 
configurations. 

The  field  variations  are  observed  to  be  considerably  different 
in  three  configurations.  The  rectangular  transducer  has  the 
largest  variations,  the  circular  next  and  the  square  on  the  di¬ 
agonal  being  the  most  uniform.  The  results  agree  with  the  ob¬ 
servations  by  Cook  and  lngenito  [10]  and  the  calculations 
along  the  axis  of  lngenito  and  Cook  [3] . 

The  variation  in  the  integrated  optical  effect  are  not  as  large 
as  the  variations  in  the  local  value  of  pressure.  Fig.  4  shows 
the  variation  in  pressure  generated  by  a  circular  transducer  of 
far  =  10  for  comparison  with  Fig.  3(a). 

Having  the  field  values  of  v  (x,  z)  allows  one  to  calculate 
photographs  of  schlieren  patterns  from  first  principles.  To 
demonstrate  how  the  transducer  geometry  and  sound  level  af¬ 
fects  the  schlieren,  we  present  a  series  of  computed  schlieren 
patterns.  The  spatial  filter  of  the  schlieren  was  chosen  to  pass 
only  the  zeroth  diffraction  order.  The  patterns  were  ealeu 
lated  by  assigning  each  point  in  relative  intensity 

W!(0  (*.*))  (13) 

where  Jo(v)  is  the  zeroth  order  Bessel  function  of  the  first 
kind.  A  logarithmic  resporse  of  the  film  was  assumed  and  a 
gray  scale  assigned.  Figs.  5,  6,  and  7  correspond  to  the  fields 
of  Fig.  3(a),  3(b),  and  3(c)  at  various  sound  levels.  In  each 
case  the  value  of  v  listed  is  the  maximum  value  found  any¬ 
where  in  the  field.  At  low  levels,  there  is  a  correspondence  be- 
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Fig.  6.  Computed  bright  field  schlicren  foi  rectangular  transducer  at 
various  sound  levels.  (a)u*2.5.  (b)c=50.  (c)u*7.J.  id)u  = 
10.0. 


ig.  5.  Computed  bright  field  schlieren  for  circular  transducer  at  vari¬ 
ous  sound  levels.  (a)ua2.5.  (b)i’“5,0.  (c)o-7.5.  (d)c=10.0. 
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Fig.  7.  Computed  blight  field  sehlieien  for  square  transducer  on  its 
diagonal,  (a)  0=2.5.  <b)o=5.0.  (c)  0=7.5.  (d)o=10.0. 


tern  shown  in  Fig.  8,  /e  have  chosen  a  spatial  filter  that  passes 
only  the  first  diffraction  orders.  Consequently,  intensity  as¬ 
signed  to  each  point  is  the  square  of  the  first-order  Bessel  func¬ 
tion  with  the  argument  o.  Other  schlieren  patterns  are  shown 
elsewhere  [11]. 


Conclusion 

Front  first  principles  it  is  possible  to  calculate  the  integrated 
optical  effect  for  any  given  transducer  configuration  in  the 
Raman-Nath  region.  The  numerical  technique  results  from 
application  for  Fourier  projection  slice  theorem  and  does  not 
require  explicit  integration.  Implementation  via  FFT  permits 
rapid  calculations  of  the  integrated  optical  effect  and  schlicrcn 
patterns. 

This  procedure  can  be  adapted  easily  to  other  configurations 
tween  the  gray  levels  and  the  field  value  of  v  in  Fig.  3.  How-  such  as  focused  transducers,  imperfect  transducers,  zone 

ever,  at  higher  sound  levels,  the  nonlinear  properties  of  (13)  plates,  and  shaded  transducers, 

tends  to  make  the  patterns  more  complicated  and  less  intcr- 
pretable.  However,  the  spreading  nature  of  the  sound  beam 
becomes  more  evident  at  the  higher  levels. 

An  increase  of  ka  also  changes  the  pattern.  From  the  de¬ 
scription  of  Osterhammcl's  experimental  arrangement  and  an 
estimation  of  the  sound  level,  wc  have  computed  the  image 
shown  in  Fig.  8  which  is  to  be  compared  with  Fig.  1 .  As  one 
can  sec,  many  of  the  geometrical  features  are  duplicated  in  the 
calculated  image.  In  fact,  the  computed  image  contains  some 
gray  areas  which  are  in  the  original  photographs  but  have  been 
lost  in  subsequent  photographic  reproductions.  For  the  pat- 


Fig.  8.  Computed  dark  field  schlieren  approximating  ihe  conditions 
of  Fig.  1  of  a  square  lransduccr  with  ka  =  45,2. 
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An  important  parameter  in  aeons to-opt it's  is  the 
so-called  Raman- Math  parameter  which  measures  the 
integrated  change  of  optical  path  produced  by  an 
ultrasonic  field.  A  st  rai  ght- forward  procedure 
has  been  developed  for  calculating  this  parameter 
for  the  entire  sound  field  for  an  arbitrary  vcloc- 
i t y  profile  and  arbitrary  configuration  of  a  baf¬ 
fled  planer  transducer,  examples  will  be  shown 
for  tiie  square  transducer.  Hie  technique  is 
limited  to  linear  acoustics  and  the  Raman-, Nath 
region  of  acousto-optics  interaction. 


figure  I  shows  the  basic  optical  arrangement  for 
visualizing  ultrasonic  fields.  The  type  of  spa¬ 
tial  filtering  peiTenned  in  plane  I',  dictates  tin¬ 
type  of  pattern  seen  in  the  image  plane  l’t.  Most 
amnion  types  of  filter  choices  are  the  following: 

1.  Pass  onlv  the  zeroth  order 

2.  block  only  the  zeroth  order 

!.  Pass  only  the  first  diffraction 
ordei 

However,  the  complexity  seen  in  the  Scbloiicn  pat¬ 
tern  depends  on  the  soundfield  and  the  integration 
of  the  optical  effect  through  this  field,  figure 
7  shows  the  complexity  that  can  occur  in  the  photo- 
giaph.  In  this  case,  the  transducer  is  rectang¬ 
ular  and  supposedly  moving  as  a  piston.  In  this 
paper,  we  demons t rate  that  these  patterns  can  he 
thenret ieal ly  pi  edict ed  by  modeling  the  configur¬ 
ation  and  movement  of  transducer,  the  opt ieal 
intcgi at  ion  over  the  resultant  sound  field  and 
the  spatial  filtering.  Ilic  mathematical  details 
of  implementing  these  calculations  are  outlined 
elsewhere'.  However,  it  should  be  noted  here  that 
the  procedure  for  calculating  the  integrated 


optical  effect  is  rapid  as  it  is  implemented  hy 
applying  the  fast  fouricr  transform  algorithms  and 
lias  no^integrat  ion  as  given  in  other  tech¬ 
niques^’  . 

Basically,  it  involves  taking  the  projec¬ 
tion  of  transducer  velocity  in  the  direction  of 
light  propagation;  taking  the  fouricr  transform, 
dividing  by  the  cosine  of  the  angles  of  plane 
wave  assigned  to  resultant  transform  variable  , 
accounting  for  propagation  with  a  phase  term  and 
taking  the  inverse  fouricr  transform.  As  these 
are  one  dimensional  transforms,  one  can  calculate 
the  field  with  close  spacing  without  excess  cost. 

In  the  calculations  shown,  we  have  applied 
simple  spatial  filtering  typically  done  .  /Ml 
arbitrary  choice  for  photographic  film  response 
has  been  included  to  yield  a  ha  lance  of  the  gray 
tones. 

figure  2  is  taken  from  a  1941  paper  hy 
Ostcrhanmel  ’.  from  the  parameter  given  and 
trial  and  error  estimation  of  the  sound  level,  we 
show  figure  A  as  an  example  of  the  calculation. 
Many  features  shown  in  Ostcrhaiiim-I  *s  photograph 
are  reproduced  in  our  computer  results. 

A  demonstration  of  the  effect  of  the  direc¬ 
tion  of  light  propagation  to  sound  field  has  pre¬ 
viously  been  published  .  It  is  duplicated  ill 
figure  4.  The  light  is  propagating  with  res¬ 
pect  to  the  oriental  inn  of  n  square  transducer 
as  illustrated,  figure  7-  shows  the  result  ol 
calculations  with  the  same  orientation  but  with 
a  different  ka. 
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IApiiuiietU.il  me  ;im  i  if  incuts  ol  acoustic  pressure  fields  near  the  focal  points  of  ultrasonic  Iransdui  eMcns 
\ stems  me  compared  to  cate nl.it ions  of  ;i  recently  developed  analytic  technique.  Ibis  meihixl  incorporaics 
i  lie  I  resile  I  approximation  u«  expand  the  Held  in  terms  of  ( iaussian  -1  agneric  functions  to  calculate  1  lie 
i  oiistie  fw  K|  lor  arbitrary  placemen  I  of  a  lens  in  llir  near  He  Id  of  a  circular  transducer.  The  predicted  field 
pallet n  sMiinuiiies  near  die  focal  region  with  a  lens  located  either  at  one  focal  distance  or.it  the  transducer 
face  arc  confirmed  by  experimental  measurement*  made  with  a  niicrosphere  probe 
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Hem  n» I v  ro|MU.*fl  results  i*f  ('av.m.iRli  and  Cook1 
«  ;l.iMi«di  a  nielliM.i  for  c.tlculal iim  Hie  acoustic  pres- 
‘-Mtc  fields  in  rttrular  ultrasonic  I  rnnschiccr— Ions  sys¬ 
tems.  I’nliko  nlhfv  fcclinirptes,  lliis  n ioI In <e|  permils 
Hie  <lt  l«  iminnt  on  of  pressure1  fields  when  a  Ions  is 
placed  at  an  aihilravv  posiliim  on  axis  in  I  hr  neaiTield 
*»f  a  rttcul.-it  |i[s(t»ii  (r.inschiccr.  These  c. denial tens 
tndicale  tli.it  sipinfic.ini  spatial  field  ehuni’os  occur  in 
Ihe  focal  region  when  I  In*  position  of  I  ho  Ions  shifts.  In 
p.o  licnl.it ,  t Ik*'  predict  spat  ial  svmntrl  t  v  of  | lie  field 
alM*iit  I  hr*  focal  nlant'  when  tin-  Ions  is  placed  mi  focal 
I'-ii’.th  'i  »m  I l*i *  l  ransclncoi ;  similarly,  wilh  llto  Ions 
l,,o.tl *  I  »i  Hit'  surface*  »»f  t h<»  transducer,  tliov  predict 
a  I'tossi  i  »*  1 1 1 .  i  >.  i  n  m  1 1 1  sfitflcd  towards  I  he  I  ranstlucor 
i  it  In* i  than  al  the  uontuol  t  ic  focal  pninl.  This  Iasi  fea- 


litre  is  a  known  fact  of  physical  oplics  and  has  been  con¬ 
firmed  experimentally  many  limes. 
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v.  in  hr*  , 


0  1b- 


UJ 

u 

?: 

< 

H- 

(/) 

a 


on! 


o  V>  *10  *"  1.5  *  2  0 

DISTANCE  f-ROM  LENS 


■’’Uorl.  p  i  linin'  l-ihth-  d  the  t  ,x  N  r x  I  Heson n  it  I ,  duo  .do i y ,  I’Kl.  J.  Cnlenhitcd  presume  liehl  lor  Ion*  iiv*  inl»-d  at  tr.in**- 

U  i.hi!v.doii,  |i.»‘  ,  inuli't  umpire**  nl  Physical  Acoustic «s  <Yn-  dui-ci  Ph  i*  tiipjioi’)  and  at  »*nn  foc.-d  d»«t.aner  lr*>ni  Ihe  ttnns- 

li-r,  f’ode  :.M»  dm  ei  f  n*o  llinvcrt.  Dimensions  tire  in  Inches 


that* 


I  AtooM  rw»r  Am  (tfMM.Nov  WHO,  0001  Wifi,  ROM  VOW  03$00  00.  t  MM0  Annul  Soc  Am  ;  Irttms  to  tin*  t  diUir 


1578 


RAD-AL  DISTANCE 


\  v,/  >, 

1 1  a-  n 


0  04  -g 


DISTANCE  FROM  LENS 


O  V 


0.B  10  12 

DISTANCE  FROM  LENS 


I 

.  * 


DISTANCE  PROM  I  FNS 


II’  I  |*  I  III*  Ml  « 1 1  \  »!:••  f  III  I  i|i|iTI  a  II- 1  I*  « 1**11  f  ill'll 

i  i  mi  miiJ.ihi  f •  *m i n* Mini c*f I  .’ll  ti  iiisdnrr  r  I  ire. 


I  f  rii'iiM’l ,ill  •niMsiirnl  values  are  here  romp.ired 
Hit  r  drill  ilrd  values  n|  Hie  t  if  If  I  in  tlir  region  near  Mu' 
Imi  ii  point  fnr  i  transducer  nf  0.5  in.  ilinni 

*pei  if  f*i»»  ii  i  lirnui  iii’V  nf  5.0  MM/.  A  Incite  Ions  with 
i  Ii"  it  !«*n»;ih  1  I  ill.  w  .is  »  .hired  rnaxinlly  linlli  .if  Mir 

1 1  .hi  d  i,i,r  f.ii'i1 11  mil  at  'Mi*-  fnp ;i|  distance  from  flip  trnn- 
fun  i  face.  I  h#  residing  pressui e  field  was  measured 
!»•.  in  111“  ,i  imernsplieri*  in  r.islerlikc  fashion  through 
ifn  |i»i*.i|  if'.inn;  iltf>  f < *r hit iia .  ( level nfii'd  I iv  Kd  wards 
.eid  -I  n/vir  Ui,  •  utili/os  Hi*  amplitude  nf  Hip  sr. dieted 
.»•  .»  up  .*•  h i < *  «»f  i|m«  lii-ld  si  length  ;if  the  ii i i «* t*« i - 
f'fi.  i«  fiM  il  mn  (In*  U  In  Mils*  an  .injimiiriil  ,  llir  mi*' in- 
pip  i*  i  p'»  ;il  mn*  if  .if  flip  Inral  pninl  »*f  Ihr  rreeiirr 
Mid  h- d Ii  (|  r  field  .il  »  IlMMf  pie.il  inn  (111  I  ill;’.  flic  I'mii  «;#» 

•  I  iti*  •  \pi  i  intent.  I  lie  I  raiisni ifter  and  h  im  system  is 
miiftl  s"  fli.il  I Im  imrinspliere  interrogates  Hip 
f»  uei|  ii  ••  if »n  *d  Mie  prrssuie  field.  The  uiiernsplie re 
In  a  t .it I ui*4  m|  appmxiinaf etv  0.00 1  in.,  suffieienf Iv 
'i'n  ill  •  e  Hi.it  even  Hip  fneused  field  presents*  essentially 
i  ti  in  f  >iin  1  i«  *  Ii  1  me  r  tin  .him  «*f  I  lie  sphere.  At  5.0 
Mil/,  Hie  iMndtiel  ««f  u*  i vermisl ant  and  diameter,  ka.  is 
■ippi  "\im,ile|v  0  '  I  he  autpl  itiide  nf  |he  siMlterrd  wave 
is  Htrirfme  if  lativrlv  independent  nf  Ihe  srallrrini: 
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I  IC.  a  M-'ioiiril  f  i  .11.1  ci>ni|xi|i*r  (-t  I'lC^^m (’  in  the  phuto 
l>. ■  i  |«*ii*liri't  1 1  in  T 1  ii«  1 1  .insniittoi  *i\is  :tt  (lie  imihl  when*  llio 
axed  i.ifp-Mif*  N  ii  in'. sninin.  Ini  llio  m  se  whore  the  lens  ip 
innuiiioil  nn  Itio  1  run sil'icc  i  T  ice. 


angles  cnriHmlorr'cl  in  llio  experiment. 

fie, mo  2  illust  rates  llio  calculated  fields  for  flic  I  wo 
1 1  ansdueor  -  Ions  geometries.  Tlio  top  diagram  illus- 
t rules  Hie  Hold  nlit .lim'd  for  llio  Ions  positioned  :il  tlio 
tv  .in  itinvi  face;  tlio  same  pattern  can  1  to  exported  for 
tlio  situation  of  i  ti.msduccr  with  a  curved  radiating 
surf. ice  llio  lower  diagram,  fur  I  hi*  Ions  positioned 
.i*  •  >iii*  focal  length.  indicates  a  sipnifieant  stiuetural 
change  noouis  undei  these  circumstances.  In  both 
i,i.  os.  Hie  aliscissa  measures  tlio  distance  from  the 
lens  and  llio  ordinate  has  boon  expanded  In  show  the 
result, ml  spatial  detail.  Minor  differences  in  llio  hit’ll 
lie". t'l'O  nl  svmmelr>.  exhibited  in  llio  lower  of  the  plots 
is  a  le  ntt  o|  tin.  is  nip'iter  graphic  einilour  routine  and 
not  t  l lit. . .iii.il v I ic  deveiopnienl;  ill  tioth  iiisl.ineos  Hie 
lire;  ane  peak  in  llio  diagrams  is'  normalised  lo  H .0  dll 
nid  the  contours  chosen  In  map  H.ll-dlT  intervals  Tlio 
i  xpei  inteiilal  and  I'.ilotilaled  values  shown  in  the  re¬ 
in. lining  figures  expand  only  the  central  region  of  those 
del  ailed  plots. 

Mi  a  ned  and  ronipulrd  pressure  profiles  for  llio 
<-.■  •  he  1 1*  Hu  ten  I  niniinlril  d  tile  1 1  .'I’l.s'ilueei  face 
an  mm  in  I  ic  a  I  lie  iiieas'ined  oniilour  slruclme 
in  tin  c  i  e  is  m  ■to. nl  ai; i eemeol  with  the  picdiefed 


values;  the  expanded  scale  of  presentation  allows  veri¬ 
fication  of  the  location  of  t he  pressure  peaks  and  nulls. 
In  Imlh  diagrams,  the  contour  nearest  the  peak  is 
-G.OdB  relalive  lo  I  lie  peak  pressure;  llio  remaining 
inlervaJs  are  -  3.0  dB  down  (i.o.,  -  0.0  dB,  -  12.0  dB, 
etc.).  There  is  an  experimental  error  of  approximately 
1.0  dB  due  to  tlio  nonlinear  response  of  the  measuring 
eleclronic  equipment. 

Figure  4  shows  llio  measured  and  computed  contours 
for  the  case  where  llio  acoustic  lens  is  placed  one  focal 
distance  from  flic  transducer  face.  The  measured  pres¬ 
sure  peak  and  nulls  are  again  in  good  agreement  will) 
Hint  of  Ihe  calculation;  the  measured  values  verify  the 
high  degree  of  symmetry  cxpeclcd.  Calculated  and  ex¬ 
perimental  pressure  profiles  along  1  lie  transducer  axis 
in  the  region  of  the  pressure  peak  are  shown  in  Fig.  5. 
The  uppor-lcf!  curve  shows  Ihe  caleulaled  pressure  for 
a  lens  moulded  on  the  source  face,  Ihe  upper  right 
shows  the  experimentally  measured  pressure.  The 
lower  curves  repeal  the  process  for  the  lens  located 
al  Ihe  focal  distance.  The  caleulaled  pressure  curve 
in  Ihe  plane  perpendicular  to  t lie  axis  al  Ihe  point  of 
maximum  axial  pressure,  a  distance  of  approximately 
0.06  in.  from  the  lens,  is  ploltcd  in  Fig.  6,  along  with 
Hie  experimentally  measured  values. 

Kxperimcnlal  evidence  indicates  that  the  method  of 
Cavanagh  and  Cook  adequalelv  predicts  tlio  pressure 
field  of  ultrasonic,  transducer-lens  system.  Calcula¬ 
tions,  performed  point  for  point  with  Hie  CJaussian- 
l.agderre  algorithm,  allow  for  quick  compulation.  The 
results  outlined  in  this  work  assume  Hie  lens  large 
compared  In  llio  width  of  Ihe  sound  Imam.  The  analytic 
formulation  lias  been  exlended  lo  predict  the  field  for 
finite  sm’  lenses.  Addition. illv.  Ihe  experiuient.il  evi¬ 
dence  presented  indicates  the  utility  oftlm  microsphero 
probe  In  spatially  map  with  high  resolution  and  preci¬ 
sion  the  fields  encounter  in  focused  ultrasound. 
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■I  Mu  i  \|>«  i  ininii .  ‘|  hi*  I i.nismiHer  .mil  lens  syMi'in  is 
it  <i»  iii.il  I  hi*  line  ii»s|ih(*n»  inlpriii(*;ili*s  llm 
iii'il  ii  imi  "l  tin*  pressm i*  fielil.  The  mic ensphere 
i  ,i  i.nhns'»l  .i|*}*io\iin;ilolv  0.001  in.,  sufficient  I  v 
mi  ill  i*  llul  even  tin1  fnensril  fiehl  piesenls  c’ssrnti;ill v 
unit  <im  hold  •»\i*r  the  ,ihm  "f  the  sphere.  At  fi.O 
HI  .  tin*  piihIui'I  nf  \v  ivefmsl.inl  :m«l  (lijmeler.  k.i.  is 
ippi  "Xim.iti'lv  O.r*.  I  he  .implilinle  <«f  llie  scntlerril  \v;ivc 
s  tin  lefiMc  ii  l.ilnelv  independent  nf  llie  senllci'int? 
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I  II.,  'i.  (‘fmiiKirlsnn  of  i*rcssiiie  imi  Iran  si  hi  nr  sivls  tlirmiuh 
l<«:i|  Vi  uh>n,  ennipntril  for  U  ns  on  tinnsihu  rr  f.n  •  hi|»i«  r 
left1,  hh  iisih'i  il  ultli  |cn«  on  fid  e  (np|MT  riuht*,  cunt  (mb -it  with 
h  ns  at  mie  fm -;i1  ili'Liiire  Hewer  hfO,  ;i n<1  on -a sined  with  lens 
;*t  nne  (oral  rilstinrc  flower  rl^htl. 
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